Time Double-Slit Interference in Tunneling Ionization 
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We show that interference phenomena plays a big role for the electron yield in ionization of atoms 
by an ultra-short laser pulse. Our theoretical study of single ionization of atoms driven by few- 
cycles pulses extends the photoelectron spectrum observed in the double-slit experiment by Lindner 
et al, Phys. Rev. Lett. 95, 040401 (2005) to a complete three-dimensional momentum picture. 
We show that different wave packets corresponding to the same single electron released at different 
times interfere, forming interference fringes in the two-dimensional momentum distributions. These 
structures reproduced by means of ab initio calculations are understood within a semiclassical model. 

PACS numbers: 32.80.Rm,32.80.Fb,03.65.Sq 



I. INTRODUCTION 

In the last years advanced laser facilities have achieved 
intensities of the order of 10 15 lV/cm 2 and pulse lengths 
of the order of 10 fs, which corresponds to few cycles of 
an electrical field of 800 nm wavelength EI- The inter- 
action process of such short and strong pulses with mat- 
ter is a topic which has attracted much interest recently. 
Many experimental (see for example 0,01) and theoreti- 
cal studies have been performed in this line. Theoretical 
calculations employ different methods: ab initio by solv- 
ing the time dependent Schrodinger equation (TDSE) 
PTIM 0- 0) or by using semiclassical approximation 
within the strong field approximation (SFA) methods 
HUES* the (Coulomb-) Volkov approximation [Til fl^| . 
or by quasiclassical approximations where the electron 
is considered classically but the possibility of tunn eling 
is incorporated within the WKB approximation |y, Il3| . 
Very recent experiments show clear signs of interference 
among different parts of the wave packet detached at dif- 
ferent times (different optical cycles of the electric field) , 
when an external linear-polarized short-laser pulse is ap- 
plied [3. Clear interference peaks in the photoelectron 
spectra are observed in this kind of experiments. The 
photoelectron spectra of He + for different angles of ejec- 
tion are calculated in jl9| . 

Another recent advance is the imaging of the momen- 
tum distribution of the ejected electron in ionization of 
rare-gas atoms by a moderately strong ultrashort laser 
pulse, in the transition regime from multiphoton to tun- 
neling ionization. In single-electron emission [3, the 
longitudinal momentum distribution (k z along the di- 
rection of the laser polarization) of photoelectrons from 
rare gases features a broad "double- hump" structure near 
threshold (E ~ 0) which surprisingly resembles the k z 
distribution for non-sequential double ionization 3] . This 
structure is due to the interplay of the Coulomb interac- 
tion and laser field on the receding trajectory of the elec- 
tron U 0, 0, 0| . A recently measured complex interfer- 
ence pattern near threshold in the two-dimensional mo- 



mentum (k z ,k p ) plane [p4lE5l| was also explained as near 
threshold Ramsauer-Townsend diffraction oscillations ■ 
A simple semiclassical analysis identified the fringes re- 
sulting from interfering paths released at different times 
but reaching the same Kepler asymptote. Very recent 
experiments show clear signs of interference among dif- 
ferent parts of the wave packet detached at different times 
(different optical cycles of the electric field), when an ex- 
ternal linear-polarized short-laser pulse is applied [l8j| . 
Non-equispaced interference peaks in the photoelectron 
spectra are observed in this kind of experiments and cal- 
culations, for example the photoelectron spectra of Hc + 
for different angles of ejection Ell- 
in this paper we generalize the study of interference in 
the energy (scalar) domain E3 to the momentum (vec- 
tor) domain for the case of a 1-2 cycle pulse. We cal- 
culate the electron yield when a hydrogen atom is sub- 
ject to a linear polarized two-cycles laser pulse with a 
sin 2 envelope function and show its similarity with a 
one-cycle pulse without envelope function. We demon- 
strate that non-equispaced peaks in the photoelectron 
spectrum stems from interference of different parts of the 
wave packets released at different times. Firstly, we solve 
the TDSE, then we make some classical considerations in 
order to intuitively understand the role of classical trajec- 
tories on interference phenomena, and lastly we present 
a semiclassical model which reproduce satisfactorily the 
calculated electron-yield distributions. 



II. THEORY AND RESULTS 

We consider a hydrogen atom interacting with a lin- 
early polarized laser field. The Hamiltonian of the system 



H 



+ V{r) + zF{t), 



(1) 



where V(r) = — 1/r is the Coulomb potential energy, 
p and r are the momentum and position of the electron, 
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respectively, and F(t) is the time dependent external field 
linearly polarized along the z direction. The laser pulse 
is chosen to be of the form 



F(t) 



f(t)sm(ujt) z (0 < t < r) 
otherwise 



(2) 



where u is the laser frequency, r the total pulse duration, 
and f(t) is the envelope function. Atomic units are used 
throughout. 



TDSE calculations 



The time-dependent Schrodinger equation is solve d by 
means of the generalized pseudo-spectral method [2Cj . 
The method combines a discretization of the radial co- 
ordinate optimized for the Coulomb singularity with 
quadrature methods to achieve stable long-time evolution 
using a split-operator method. It allows for an accurate 
description of both the unbound as well as the bound 
parts of the wave function \ip(t)). The process of detect- 
ing an electron of momentum k can then be viewed as a 
projection of the wave function onto the Coulomb wave 
functions 0,0, Therefore, after the laser pulse is 
turned off, the wave packet is projected onto outgoing 
Coulomb wave functions which gives the asymptotic mo- 
mentum distributions 



dP 
dk 



1 



i 



J8i(k) ^2l + lPi(coa9 k ) (k,l |^(r)) 



(3) 

where 5i(k) is the momentum dependent Coulomb phase 
shift, 9k is the angle between k and the polarization di- 
rection of the laser field, z, Pi is the Legendre polynomial 
of degree I, and \k,l) is the eigenstate of the free atomic 
Hamiltonian with positive eigenenergy E — k 2 /2 and 
orbital quantum number I. The Coulomb projection is 
needed for the cases that a physical magnitude is not a 
constant of motion of the free evolution (once the external 
field is turned off), i.e., the components of the kinetic mo- 
mentum. In turn, this is not the case for the photoelec- 
tron spectrum since the energy is a constant of motion of 
the free evolution. Cylindrical symmetry makes the dy- 
namics a two dimensional problem, where the projection 
of the angular momentum on the polarization direction 
of the laser is a constant of motion (the magnetic quan- 
tum number m is unaffected during the time evolution). 
As the initial state of the system we consider the ground 
state of the Hydrogen atom, (m = 0). 

Firstly, a driving field of frequency to = 0.05 with enve- 
lope function f(t) = F sin 2 (7rf/r), peak field F„ = 0.075 
corresponding to an intensity of 2x 10 14 W/cm 2 , and time 
duration r = 251 (which involves two complete optical 
cycles) is used. The inset of Fig. 1 (a) shows the pulse 
shape. The photoelectron spectrum, shown in Fig. 1 
(a), has a non-equally spaced peaked structure, whose 
separation between two consecutive peaks increases with 



energy. Recent simulations 

.19] have also shown the ex- 
istence of these peaks in the density probabilities at dif- 
ferent angles of ejection but for a much stronger field 
(10 16 Wcm~ 2 ) applied to He + atoms calculations. 
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FIG. 1: (a) Photoelectron spectra for a two-cycle 
electric field of frequency to = 0.05 a.u., peak field 
Fq = 0.075 a.u. and duration r = 4ir/ui = 251. Inset: 
Electric field as a function of time, (b) Doubly differ- 
ential electron momentum distribution in cylindrical 
coordinates {k z , k p ) in logaritmic scale for the same laser 
pulse. 

In order to develop a more complete analysis we show 
in Fig. 1 (b) the double differential momentum dis- 
tribution as a function of the final longitudinal — k z — 

and transversal — k n = » kl + k% — momentum of the 



electron 



dk dk z 



Different characteristics of the two- 



dimensional distribution are worth to be highlighted: 
(i) about the 95% of the distribution lies in the region 
of positive longitudinal momentum (k z > 0), (ii) more 
specifically the distribution is constrained to the region 
—0.5 < k z < 2.5 and k p < 0.6, and the most important 
,(iii) a pattern with strip shape. It can also be observed 
some small radial structure for positive k z near thresh- 
old. This kind of structure has recently be understood as 
generalized Ramsauer-Townsend type diffraction oscilla- 
tions resulting from interfering paths released at different 
times 0- 

In the following we consider a pulse with constant en- 
velope function f(t) = F within [0,t] and outside for 
two different values for the pulse duration: r = tt/lu = 63 
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(half-cycle pulse) and r = 2tt/lo = 126 (full-cycle pulse). 
In the case of the half-cycle pulse, the field, obviously, 
does not correspond to a traveling but to a standing 
wave [see Fig. 2 (a)]. We can observe the similarity 
of the shape of the one-cycle pulse of Fig. 2 (a), with 
the two-cycle pulse with envelope function in the inset of 
Fig. 1, where we can consider only one effective cycle (at 
the center) while the other is rather weak. As the pulse 
shapes of these two pulses are similar so should it also be 
the electron yield distributions. 
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FIG. 2: (a) Electric field of frequency u = 0.05, 
duration r = 2it/lu = 126 a.u. and amplitude 
F = 0.075 a.u. as a function of time (1 cycle and 
square envelope function), (b) Kinetic momentum for 
different classical trajectories: electron detached at time 
to = 0, 7r/2w, tt/uj, 37r/2w, and 2ti/lo, as indicated. The 
momentum of the trajectories to — tt/2uj and 3tt/2uj 
coincide in the time interval t > 3tt/2lu. 

In Fig. 3 (a) the photoelectron spectra for one-cycle 
pulse must be compared to the one of Fig. 1. In Fig. 3 
(a) we showed the photoelectron spectra for a half-cycle 
and a full-cycle pulses are shown. The total ionization 
probability, which can be calculated as the integral of the 
energy distribution, is higher for a full cycle pulse than 
for a half cycle pulse, as expected. After half a cycle 
the ionization probability is 0.06, while after a full cycle 
pulse it is the double, 0.12. This shows that the system 
is in the linear ionization regime very far away from the 
saturation limit (depletion of the ground state). We will 
make use of this fact below, in subsection C. Both energy- 
distributions feature a peak structure with bigger separa- 
tion between consecutive peaks as the energy increases. 
In turn, the photoelectron spectrum of a half-cycle pulse 
shows no oscillating structure. In Fig. 3 (b) the double- 
differential momentum distribution for a half-cycle shows 



a complete smooth distribution in the negative longitudi- 
nal momentum -k z - region. The pulse is not only respon- 
sible for the ionization of the atom but also for the motion 
of the wave packet towards the negative longitudinal mo- 
mentum region. In turn, when a second half-cycle is in- 
cluded [Fig. 3 (d)] the whole wave packet moves towards 
the positive region and a clear fringe structure appears. 
This fact leads us to think that the peak structure of 
the one-cycle pulse spectrum emerges as a consequence 
of an interference effect. When considering a full cycle 
pulse, the wave packet created (released) during the first 
half cycle is pulled to the negative longitudinal momen- 
tum region after the first half cycle, but then is pulled 
towards the positive longitudinal momentum region by 
the second half-cycle. We observe that while the wave 
packet is spread about 3 a.u. in the positive longitudinal 
momentum, about only 0.5 a.u. (17%) in the transversal 
momentum. Additionally to the ionization in the first 
half cycle, another part of the wave packet appears in 
the continuum during the second half cycle which is also 
pulled towards the positive longitudinal momentum re- 
gion. These two wave packets interfere, which can be 
observed in the double differential momentum distribu- 
tion as clear interference fringes [Fig 3 (c)] and in the 
photoelectron spectrum as non-equispaced peaks [Fig. 3 
(a)]. The spacing between the peaks increases with en- 
ergy. 
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FIG. 3: (a) Ab initio photoelectron spectra for the 
case of the electric field in Fig. 2 (a): after a half-cycle 
pulse (r = tt/uj) and one-cycle pulse (r ~ 2it/lu). (b) 
Ab initio doubly differential electron momentum distri- 
bution after the half-cycle pulse (one slit), (c) Ab initio 
doubly differential electron momentum distribution after 
the one-cycle pulse (two slits). 
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These interference processes can be reproduce with a 
simple semiclassical model of two electronic wave packets 
escaping from the nucleus near the extremes of the elec- 
tric field and interfering in momentum space. We will 
see in subsection C that the mentioned interference pro- 
cess can be understood as the interference in momentum 
space of two electron classical trajectories escaping from 
the nucleus, one at each half cycle near the extremes of 
the electric field. We will give a clear analysis of this 
interference effect below by means of a semiclassical ap- 
proach, but before a few classical considerations should 
be stated. 



B. Classical considerations 

In this subsection we want to intuitively understand 
the physics of interference described in last subsection. 
We will see that it is sufficient to use a one-dimensional 
picture to glimpse the interference of the ionized atom. 
By using the classical equation of motion and considering 
the electron to detach from the nucleus at time to with 
zero velocity (Simple Man's Model [23J), it acquires the 
following (longitudinal) momentum 



same final momentum p\ — P2, if their respective release 
times t\ and ti accomplish with the relation 

w (h + ta) = 2vr. (5) 

This is exactly the condition for interference of two clas- 
sical trajectories in the semiclassical model, as we will 
see in the next subsection. 

In Fig. 2 we have considered only "direct electrons", 
that is those which do not suffer any process of rescatter- 
ing. In this one-dimensional model a rescattering process 
means an inversion of ir in the scattering angle (change of 
sign in the momentum). Thus, the electrons would travel 
to the negative position direction, —z, instead of being 
expelled towards the positive direction after a complete 
cycle. As there are no other trajectories going in the — z 
direction, so there is no possibility of interference. There- 
fore, next these "indirect electrons" can be neglected for 
a one-cycle pulse and will be excluded from our semi- 
classical model of interference. This is not the case for 
a longer (more than one cycle) pulse where rescattering 
can occur in different cycles with the ensuing production 
of interference arising from indirect electron trajectories. 

C. Semiclassical model 



Fq 

p(t) = — (cos ujt — coswio) 

UJ 



(4) 



In Eq. we make use of the strong field approxima- 
tion (SFA), thus the effect of the atomic potential on the 
detached electron is neglected. The momentum p corre- 
sponds to the longitudinal momentum k z in the three- 
dimensional case. 

Most probably the atom is ionized at the extremes of 
the one-cycle electric field, i.e., ujt a = tt/2 and 37r/2, thus 
Eq. 10} leads to p(t) = (F /oj) cos u>t. Therefore, after a 
half-cycle pulse (r = tt/lu) the momentum of the electron 
is p(t) = —Fq/oj, and after a complete-cycle pulse (r = 
2-k/u), p(t) = F /lj [see Fig. 2 (b)]. In Fig. 3 (b) and 
(c) we can observe that the center of mass of the wave 
packet follows the classical quiver motion with amplitude 
Fo/uj, i.e., right after the first half cycle it is negative, 
(k z ) ~ —Fq/lu — —1.5, and right after the complete cycle 
it is positive, (k z ) ~ Fq/uj = 1.5. The limit cases are 
trajectories of the electron released right at the beginning 
(or end) of the pulse to = (or 2tt/lu) and right after 
the first half cycle (to = tt/uj). From Fig. 2 (b) [and 
from Eq. (0}] it can be easily seen that the momentum 
is constrained to values —2Fo/oj < p < 0, right after 
a half-cycle pulse, and < p < 2Fo/uj, right after a 
one-cycle pulse. In Fig. 3 (b) and (c) we observe that 
the quantum wave packets follow this classical constrains 
quite accurately. 

In Fig. 2 (b) we see an example of two different tra- 
jectories released at times to = 7r/2w and 37r/2w having 
the same final momentum. It is straightforward from Eq. 
(0} to deduce that every pair of trajectories will have the 



Throughout this subsection we solve the TDSE under 
a set of approximations and will arrive to a simple ana- 
lytical semiclassical solution 24] . These approximations, 
some of them already mentioned, are to be introduced 
throughout this subsection. We consider that there is 
only one atomic bound state (|0)) neglecting the bound 
excited states. In order to solve the Schrodinger equa- 
tion with Hamiltonian given by Eq. Q we assume the 
following ansatz 



a(t) |0) + / dvb(v,t) \v) 



(6) 



where \v) is an eigenstate of the continuous spectrum 
with velocity v, and I p is the atomic ionization potential 
(Ip = 0.5 for the hydrogen case). When we include Eq. 
© into the Schrodinger equation we arrive to the follow- 
ing expression for the amplitude of the continuum states: 



b(k,t) = -j2 



1=1 



2mF( J t l SP ) 



1/2 



\k + A(V SP )\ 
where S accounts for the Volkov action, 

~(k + A(t")) 2 



d*(k+A(t l SP ))e lS ^sp' t \ 



S(t',t) = 



dr 



(7) 



(8) 



d*(v) — (v\z\Q) is the dipole element of the bound- 
continuum transition, A(t) = — J Q * dt' F(t') is the vec- 
tor potential of the laser field divided by the speed 
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of light, and k is the canonical momentum defined as 
k = v(t) — A(t). To arrive to Eq. J7J) we have made 
use of (v \v') — S(v — v 1 ), and the following approxima- 
tions: (i) we neglect the depletion of the ground state 
(a(t) — 1), (ii) continuum-continuum transitions are of 
the form (v\ z {v 1 ) = iV„c)(z7— v"), and (iii) in order to get 
an analytical solution we apply the saddle-point approx- 
imation method [2S|. t l SP (i = 1,2) are the solutions of 
the stationary phase action 



direct consequence of the assumption that no depletion 
of the ground state exists. 



dS(t' = t SP ,t) 



k + A{t SP ) 



dt> 2 
The solutions of Eq. @ are complex: 



I P = 0. 



(9) 



1 sp — 



L SP — 



1 



2tt 
t 

OJ 



1 



(k z T i J' 



2I„ 



k 2 p 



u> 



(10) 



l 

sp- 



Therefore, the sum in Eq. J2J must be understood as a 
two-term sum over the two different times t\ P and t 2 SP 
of Eq. UlUf l. These two terms correspond to two different 
classical trajectories. The two ionization times in Eq. 
(I10JI accomplishes with the interference condition of Eq. 

Finally, the photoelectron spectrum can be written as 
g = 2n f^d cos 9 k V2E b(k , t = 2tt / u) 

ergy is a constant of motion of the free evolution (after 
the electric pulse is turned off), thus the energy distribu- 
tion is invariant when taking the asymptotic limit t — > oo. 
In this case, the kinetic momentum is not a priori invari- 
ant, but within the strong field approximation, once the 
electron is detached, it is not affected any longer by the 
core potential and, in consequence, the momentum dis- 
tribution does not change when taking the limit t — > oo. 



since the en- 



The momentum distribution ^£ = b(k.t = 2tt/uj) 
be written as 
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FIG. 4: (a) Semiclassical photoelectron spectrum 
corresponding to the field in Fig. 2 (a) . (b) Semiclassical 
doubly differential electron momentum distribution after 
the pulse. 



6(fc,i = 27r/w) = B(k) cos 2 [AS (k)} 



(11) 



where the phase AS(k) = S(t 2 SP ,t) - S(t^ P , t) is respon- 
sible for the interference process, and the factor B(k) is 
the ionization probability at a time t l SP and states only 
for the modulation of the distribution of the unbound 
electrons. By making use of the simplified version of the 
saddle point method, where the ionization times of Eq. 
(|10(l are considered real (Simple Man's Model), the ion- 
ization probability can be written as [19| 



B{k) 



2(2I p + k 2 p )\F(t(k))[ 



■ exp 



2(27 p + k 2 ) 



3\F(t(k))\ 



(12) 

The reader can easily check that B(k) is invariant under 
the interchange tg P < — > t 2 SP in Eq. i|12|) . which is a 



The results of the semiclassical model [Eq. (|llf) ] are 
shown in Fig. 4. We observe a good qualitative agree- 
ment with the ab initio results of Fig. 3. The photoelec- 
tron spectrum calculated with the semiclassical model in 
Fig. 4 (a) reproduces the oscillating structure of the ab 
initio calculations in Fig. 3 (a), but the semiclassical 
model fails near threshold and in the high energy region 
(near classical boundaries). The same observation can be 
performed for the doubly differential momentum distri- 
butions in Fig. 4 (b) for the semiclassical model and Fig. 
3 (c) for ab initio calculations. The interference pattern 
happens to be more complicated for ab initio calculations 
even though the semiclassical model reproduces the over- 
all interference fringes. While the semiclassical distribu- 
tions sticks (by construction) to the classical boundaries, 
the ab initio distribution invades a little bit the classi- 
cally forbidden zone. 
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FIG. 5: (a) Photoelectron spectra by solving the 
TDSE and using the semiclassical (sc) model corre- 
sponding to the field in Fig. 2 (a), (b) Longitudinal 
momentum distributions by the two methods, (c) Zoom 
of Fig. (b) to show the agreement in separation be- 
tween two consecutive peaks, even though there is a shift. 

In order to perform a more quantitative comparison we 
show in Fig. 5 the results of both ab initio and semiclas- 
sical calculations. For this purpose, we have removed the 
excited bound states from the hydrogen spectrum in the 
quantum calculation (which is present in the semiclassi- 
cal model). Therefore, the interference peaks in the pho- 
toelectron spectrum and the momentum distribution are 
more enhanced since they are not masked by the ioniza- 
tion coming from other bound states. Nevertheless, the 
distributions of Fig. 5 are very similar to the one with 
the full spectrum (Fig. 3). In Fig. 5 (a) we observed that 
for a long range of energies the agreement in the position 
of the peaks of the semiclassical energy spectrum and the 
TDSE one is very good. According to the classical me- 
chanics [Eq. J3J] the energy spectrum is constrained to 



E < 2 (-^j) = 4.5 a.u. and, as mentioned before, we 
observe in Fig 5 (a) that the agreement is not so good 
near the extremes of the energy range, were the semiclas- 
sical spectrum decreases more abruptly. This feature is 
also observe for the longitudinal momentum distribution 
4rr- in Fig. 5 (b). In most of the range the agreement 
is good as it can be seen in a zoom of the momentum 
distribution [Fig. 5 (c)], but it is not very good near 
the borders of the classical domain < k z < = 3 
a.u., as expected according to the semiclassical theory. 
In both energy and momentum distribution the envelope 
function of Eq. (|T2"|> does not reproduce accurately the 
quantum results. Nevertheless, the interference pattern 
of ab initio calculations can be reproduced by the semi- 
classical model given by the second factor of Eq. Qlip. 
In figures 5 (a) and (c) the distance between two consec- 
utive peaks is very similar in the two models but there 
is a shift in the position of the peaks of one model with 
respect to the other. This shift could be due to the effect 
of the coulomb potential on the ejected electron, which 
is present in the ab initio calculations but is neglected in 
the semiclassical model (SFA). 



III. CONCLUSIONS 

In this article we have shown theoretical studies on the 
interference effects observed in the electron distributions 
of ionized hydrogen atoms subject to a linearly polarized 
short laser pulse. We have extended previous analysis 
in the energy domain to a full three-dimensional one in 
the momentum domain. The two-dimensional electron 
momentum distribution after a full cycle pulse evidences 
interference fringes. We have recognize the peak struc- 
ture in the photoelectron spectrum and the longitudinal 
momentum distribution as the interference phenomenon 
between the wave packets released at the first and sec- 
ond half cycle, functioning each half cycle as an indepen- 
dent slit. In the understanding of the interference phe- 
nomenon we have made use of a simplified semiclassical 
model which has the advantage of being analytical. The 
semiclassical model reproduces quite well the interference 
patterns in the spectra of ejected electrons. 
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